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Abstract 

We investigate the problem of transition from galactic cosmic rays to extragalactic 
ultra high energy cosmic rays. Using the model for extragalactic ultra high energy 
cosmic rays and observed all-particle cosmic ray spectrum, we calculate the galactic 
spectrum of iron nuclei in the energy range 10^ — 10^ GeV. The flux and spectrum 
predicted at lower energies agree well with the K AS CADE data. The transition 
from galactic to extragalactic cosmic rays is distinctly seen in spectra of protons 
and iron nuclei, when they are measured separately. The shape of the predicted iron 
spectrum agrees with the Hall diffusion. 
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1 Introduction 



There are now convincing evidences that ultra high energy cosmic rays (UHECR) 
in the energy range 1 x 10^^ — 8 x 10^^ eV are extragalactic protons. These 
evidences include: (i) Measurement by HiRes detector [1] of x^ax in EAS lon- 
gitudinal development favours protons as primaries at E > 1 x 10^^ eV (see 
Fig. 1), (ii) The energy spectra measured by Akeno-AGASA [2], Fly's Eye [3], 
HiRes [4] and Yakutsk [5] detectors clearly show the dip [6] (see Fig. 2) , which 
is a signature of interaction of extragalactic protons with CMB radiation, (iii) 
The beginning of the GZK cutoff up to E' ~ 8 x 10^^ eV is seen in all data, 
including that of AGASA (see Fig. 2). 
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Fig. 1. The HiRes data [1] on mass composition (preliminary) , The measured Xjjjax 
at E > Ix 10^^ eV are in a good agreement with QGSJet-Corsika prediction for 
protons. 

It is interesting to note that the excess, detected by AGASA [7] at ~ 10^^ eV 
from directions to the galactic sources, Galactic Genter and Gygnus, may be 
naturally interpreted as evidence of diffuse extragalactic flux at £^ > 1 x 
10^^ eV. Indeed, if transition from galactic to extragalactic diffuse flux occurs 
due to failure of magnetic confinement in the Galaxy, the unconfined flux from 
the galactic sources should become visible above energy of the transition. 

On the other hand, not all experimental data agree with pure proton composi- 
tion at > 1 X 10^^ cV. While data of HiRes-MIA[8] and Yakutsk [9] support 
such composition, the data of other detectors, such as Fly's Eye [10], Haverah 
Park [11] and Akeno [12] favour the mixed composition at £^ > 1 x 10^^ eV. 
The indirect confirmation of the proton composition at > 1 x 10^^ eV comes 
from the KASCADE data [13] , which show disappearance of the heavy nuclei 
from the CR flux at much lower energies. 

UHEGR have a problem with the highest energy events. First of all, 11 AGASA 
events at £■ > 1 x 10^° eV comprise signiflcant excess (not seen by HiRes) 
over predicted flux (see Fig. 2). This excess may be interpreted as the new 
component of UHEGR, probably due to one of the top-down scenarios (see 
[14]). Even if to exclude the AGASA data from analysis, the problem with 
particles at £^ > 1 x lO^'' eV remains. There are three events with E > 
1 X 10^° eV observed by other detectors, namely, one Fly's Eye event with 
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Fig. 2. Comparison of calculated spectra for non-evolutionary model (see the 
text) with observational data. There is a good agreement of a dip centered at 
£' ~ 1 X 10^^ eV with all data. The dip is produced due to pair-production 
P+7CMB p+e'^ + e~ on CMB radiation. AG AS A excess needs for its explanation 
another component of UHECR (shown by dashed curve), which can be due to one 
of the top-down scenarios [14]- 

E ^3x10^° eV, one HiRes event with E Ri 1.8 x 10^° eV and one Yakutsk 
event with ^ 1.0 x 10^° eV. The attenuation length of protons with energies 
(2 - 3) X 10^0 eV is only 20 - 30 Mpc. With correlations of UHECR directions 
with AGN (BL Lacs) observed at (4 - 8) x lO^^^ eV [15], the AGN must be 
observed in the direction of these particles. If to exclude the correlations with 
AGN from analysis, the propagation in the strong magnetic fields becomes 
feasible [16]. Nevertheless, also in this case the lack of the nearby sources in 
the direction of highest energy events (e. g. at £" ~ 3 x 10^° eV) remains 
a problem for reasonable local field ~ 1 nG. Indeed, the defiection angle 
^ ~ ^att/f^H = 3.7°ifnG IS Small, and the source should be seen within this 
angle. 

Another problem with UHECR exists at its low energy edge. The transition 
from galactic to extragalactic cosmic rays may occur at position of the second 
knee. The energy of the second knee varies from 4 x 10^'' eV to 8 x 10^^ eV 
in different experiments (Fly's Eye - 4 x 10^''' eV [3], Akeno - 6 x 10^''' eV 
[17], HiRes - 7 x 10^^ eV [4] and Yakutsk - 8 x 10^^ eV [5]). This energy 
is close to ~ 1 X 10^^ eV, where according to the HiRes data [8] protons 
start to dominate the fiux. On the other hand, KASCADE data [13] show 
disappearance of galactic cosmic rays at much lower energy. In Fig. 3 one 
can see that positions of the proton, helium, carbon and iron knees can be 
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Fig. 3. The calculated spectra of extragalactic protons (thick dashed curve labeled 
"extr. p") and of galactic iron nuclei (dash-dotted thick curve labeled "Fe model") 
in comparison with the Akeno and KASCADE data. The values of jg and used 
in the calculations are indicated in each panel. The all-particle spectrum of Akeno is 
shown by filled triangles and by solid curve. The sum of the curves "extr. p" and "Fe 
model" is exactly equal to all-particle spectrum of Akeno, given by curve "AKENO". 
The data of KASCADE for different nuclei are shown as : protons - by filled circles, 
helium - by diamonds, carbon - by stars, and iron - by inverted triangles. The 
thin full curves smoothly connect the data points. The arrows labeled by p. He, C 
and Fe show the positions of corresponding knees, calculated as Ez = ZEp, with 
Ep = 2.5 X 10® GeV. One can notice the satisfactory agreement between calculated 
and observed positions of the knees. 

tentatively accepted as 2.5 x 10^^ eV, 5.5 x 10^^ eV, 1.8 x 10^'' eV, and 
7.0 X 10^® eV, respectively (see arrows in Fig. 3). 

How the gap between the knee for iron nuclei and the beginning of the extra- 
galactic component is filled? Why the transition reveals itself in the form of 
hardly noticeable feature in all-particle spectrum? 
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2 The model 



We will construct here the phenomenological model, which predicts the spec- 
trum of iron nuclei in the energy range 8 x 10^ — 1 x 10^ GeV. Following 
Ref. [6], wc will calculate the spectrum of extragalactic protons in the energy 
range 10*^ — 10^^ GeV. Subtracting this proton spectrum from the observed 
all-particle spectrum (we shall use Akeno-AGASA data), we obtain residual 
spectrum, which, inspired by the KASCADE data, we assume to be comprised 
by iron nuclei. This spectrum will be compared with that measured by KAS- 
CADE. The justification of this model consists in confirmation of the dip in the 
calculated proton spectrum by observations. One of the results is description 
of transition from galactic to extragalactic cosmic rays. The predicted physical 
quantity is the flux of iron nuclei and their spectrum at 10^ — 10^ GeV. 



2. 1 Spectrum of extragalactic protons 



We will calculate the extragalactic proton spectrum, following Ref. [6], in the 
model with the following assumptions. 

We assume the generation spectrum of a source 

Lp{z) 
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with 
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at Ea < Ec 



E-\Eg/E^)-^o at Eg>E^ 



(2) 



The diffuse spectrum is calculated as 

ME) = ^ ^^l^ J dtqg.. (EgiE, z), E) ^ , (3) 



where Cq — TigLp is emissivity (with rip and Lp being the comoving density of 
the sources and luminosity, respectively), Eg{E,z) is generation energy of a 
proton at epoch z, dEg/dE is given in Refs. [18,19] as 

'''' [l + zg)exp{ — I dz^=4=^=—[—^) 

"-^ / E'={l+z)Eg{E,z) 



dE ' \HJ^ ^Q^(l + ^)3 + fiA 
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where bo{E) is the proton energy loss dE/dt ed, z — 0, and dt is given as 



with Hq, and I^a being the Hubble constant, relative cosmological density 
of matter and that of vacuum energy, respectively. 

In calculation of the flux according to Eq. (3) we use the non-evolutionary 
model, taking Cq not dependent on redshift z. Since we are interested in 
energies E >1 x 10^ GeV the maximum redshift appears in Eq. (3) automat- 
ically due to sharp increase of Eg{z) at large z. The control calculations with 
-^max = 4 give the the identical results. Notice that in evolutionary models the 
calculated flux is sensitive to Zuiax- 

The calculated spectrum is displayed in Fig. 2 with Ec = I x 10^ GeV and 
7g = 2.7 in comparison with observational data. In the further calculations we 
fix 7(, = 2.7 as the best fit to the observational data displayed in Fig. 2 and 
consider Ec as a free parameter. 

2.2 Galactic iron spectrum and transition from galactic to extragalactic cos- 



In Fig. 4 we show the iron galactic spectrum found as subtraction of the 
extragalactic proton spectrum calculated above, from all-particle spectrum 
(Akeno-AGASA). The transition from galactic to extragalactic cosmic rays 
occurs at crossing of proton and iron spectra, which is very noticeable and 
occurs at energies 1.5 x 10*^ GeV, 3.6 x 10* GeV and 6.2 x 10® GeV for the 
curves 1-1', 2-2' and 3 - 3', respectively. In Fig. 5 the fraction of iron nuclei 
in the total fiux is shown as function of energy. Notice, that in our model in 
energy range (1 — 10) x 10® GeV the total flux is comprised only by protons 
and iron. The Haverah Park data [11] confirm that it can be the case. The 
analysis of the Haverah Park data has been performed for energy-independent 
composition, however comparison of the obtained composition (34% of protons 
and 66% of iron) with our average value favours large Ec- 

The fiux of iron nuclei calculated here fit well the KASCADE data at energy 
E ~ 10® GeV and below (see Fig. 3). 

Why there is no pronounced feature in the total spectrum which corresponds to 
transition from galactic to extragalactic component? Such feature is inevitably 
faint when both components have approximately equal spectrum exponents. 
This is our case: the Akeno spectrum up to 7 x 10^ GeV is characterised 



dt = 



dz 



(5) 
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Fig. 4. Calculated spectrum of extragalactic protons (curves 1, 2, 2>) and of galactic 
iron spectra (curves l',2',3'J compared with all-particle spectrum from Akeno and 
AG AS A experiments. The galactic iron spectrum is obtained by subtraction of the 
calculated proton spectrum from the all-particle spectrum. The pairs of curves 1 
and 1' , 2 and 2', 3 and 3' correspond to equal to 3 x 10^ GeV, 1 x 10^ GeV, 
and 2 x 10^ GeV, respectively. The intersections of the curves 1 — 1', 2 — 2' and 
3 — 3' give the transition from galactic (iron) to extragalactic (protons) components, 
which occurs at 1.5 X 10^ GeV, 3.6 x 10^ GeV and 6.2 x 10^ GeV, respectively. 
The KASCADE data are shown by filled squares for all-particle fluxes and by open 
circles - for iron nuclei fluxes. 

by 7 = 3.0, while at higher energies - by 7 = 3.25. However, one can see 
from Fig. 3 that intersection of extragalactic proton spectrum ("extr. p") and 
galactic iron spectrum (curve "Fe-model" ) is characterised by quite prominent 
feature. Therefore, measuring the spectra of iron and protons separately, one 
can observe the transition quite distinctively. 



2.3 Visibility of galactic sources at E >1 x 10^ GeV 

We argued above that at energy E > Ec ^ 1 x 10^ GeV the extragalactic flux 
dominates. However, as the AGASA observations show [7] at these energies 
the excess of events in the direction of some galactic sources is observed. We 
shall show in this section that visibility of galactic sources at E > Ec is a. 
natural prediction of our model. 
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Fig. 5. Fraction of iron nuclei in the total flux as function of energy. 

If galactic sources accelerate particles to energies higher than 1 x 10^ GeV 
their "direct" flux can be seen, while the produced diffuse flux should be 
small because of short confinement time in the galaxy. If generation spectrum 
is dominated by protons, the "direct" flux must be seen as the protons, while 
the diffuse galactic flux is presented by the heaviest nuclei. 

We shall argue here that galactic point sources should be observed as the 
extensive sources, as it occurs in the AG ASA observations [7] . 

Propagation of the protons from the galactic sources takes place in quasi- 
diffusive regime with the large diffusion coefficient. The protons arrive at an 
observer along a trajectory in the regular galactic magnetic field, while the 
small-angle scattering in the random magnetic field provides the angular dis- 
tribution of protons in respect with the main direction along the trajectory. 
Because of the large proton energy, the scattering occurs mainly on the basic 
scale of the turbulent (random) magnetic field Ic ~ 100 pc. The scattering 
angle 9 corresponds to {9'^) = 11/ r\, where is the Larmor radius in the 
basic magnetic field. Using {B\) — (2/3) (S^) for the field perpendicular to a 
particle trajectory one obtains 

{0') = le'm)lE\ (6) 
where e is electric charge of the proton. 

After n scattering in the random magnetic field the angle with respect to the 
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direction of the regular trajectory becomes [21,22] (0^) 



Using n = r/lc, where r is a distance to the source along the regular trajectory, 
one finally obtains 



similar to formulae used in Ref. [23,24]. 

The distribution of arriving particles over the angles is the normal Gaussian, 



with dispersion a = For a distance r — 10 kpc, B — 1 fj,G and 

E = 1 X 10^ GeV one obtains a — 25° in a good agreement with the AGASA 
observations. 

At larger energies and smaller distances the sources are seen as the point-hke 
ones. 



3 Discussion and conclusions 

One of the assumptions of our model is that at high energies E > Ix 10^ GeV 
only iron nuclei as galactic component survive. This assumption is inspired 
first of all by the KASCADE data (see Fig. 3) which show how hght nuclei 
gradually disappear from the spectrum with increasing the energy. This result 
is in the accordance with other experiments, e.g. [20], which also demonstrate 
that the mean atomic weight {A) (and hence {Z)) increases with energy. 

The good agreement of calculated iron spectrum with that measured by KAS- 
CADE at £" < 10^ GeV might be delusive. The published spectra are still 
preliminary and they are slightly different in different publications. In effect, 
it is more proper to speak only about rough agreement of calculated spectrum 
with one measured by KASCADE. On the other hand, we have a free param- 
eter Ec with help of which we can fit our spectrum to the KASCADE data. 
Such fit would restrict our predictions for the Fe/p ratio at higher energies 
(see Fig. 5). 

In principle there are three classes of models explaining the knee (for a review 
and references see [25]): 

(i) Rigidity-dependent exit of cosmic rays from the Galaxy, e. g. [26,27,28]. 

(ii) Rigidity-dependent acceleration [29,30,31]. 

(iii) Interaction-produced knee [32,33]. 
For more references see [25]. 




(7) 
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As far as the spectra are concerned, the first two classes can result in identical 
conclusions. For example the spectra of nuclei in Ref. [29] with acceleration 
in presupernova winds have the same rigidity dependence as in the models of 
class (i). We shall restrict ourselves in further discussion by models of class 

(i) - 

The data of KASCADE (Fig. 3) are compatible with rigidity dependent knees 
for different nuclei, as it must be in models (i) and as it might be in models 

(ii) . 

The positions of the proton, helium, carbon and iron knees in the KASCADE 
data can be tentatively accepted as 2.5 x 10^ GeV, 5.5 x 10^ GeV, 1.8 x 
10^ GeV, and 7.0 x 10^ GeV, respectively (see arrows in Fig. 3). The accuracy 
of the knee positions is not good enough and is different for different knees. 
While position of the proton knee, Ep f« 2.5 x 10^ GeV, agrees with many other 
measurements (sec [25] for a review), the position (and even existence) of the 
iron knee is uncertain. The energies given above should be considered not more 
than indication. The future Kascade- Grande [34] data will clarify the situation 
with the iron knee. However, it is interesting to note that the positions of the 
nuclei knees, taken above from the KASCADE data, coincide well with simple 
rigidity model of particle propagation in the Galaxy. In this model Ez = ZEp, 
where Z is a charge of a nucleus. Taking Ep = 2.5 x 10^ GeV, one obtains for 
helium, carbon and iron knees -Ehc = 5.0 x 10® GeV, Ec = 1-5 x 10^ GeV and 
Epe = 6.5 X 10^ GeV, respectively, in a good agreement with the KASCADE 
data cited above (see Fig. 3) . 

In the rigidity-dependent models (i) one can predict the spectrum shape below 
and above the knee. Below the knee, in the case of the Kolmogorov spectrum 
of random magnetic fields the diffusion coefficient D[E) oc E^^'^, and the 
generation index for galactic cosmic rays is 7^ = 2.7 — 0.33 ^ 2.35. If above 
the knee D{E) oc E'', the spectral index 7 = 7c, + /c. At extremely high 
energy, when the Larmor radius becomes much larger than the basic scale 
of magnetic field coherence Ic, the diffusion reaches the asymptotic regime 
with D[E) oc E"^. In this case 7 = 4.35. The critical energy is E^r = ZcBqIc ~ 
2.4 X 10^ GeV for iron (Z=26), where Ic — 100 pc and the regular magnetic 
field on the basic scale is Bq ~ IfiG. Therefore, the asymptotic regime with 
7 = 4.35 is valid at ^ 2.4 x 10^ GeV and in the energy region of interest 
the intermediate regime with non-power law spectrum may occur. 

Below the knee the longitudinal diffusion (along the regular field Bq) domi- 
nates. It is provided by condition of strong magnetizing. At higher energy the 
role of transverse diffusion might be important. An example of it can be given 
by the Hall diffusion [26,22], associated with the drift of particles across the 
regular magnetic field. In this case above the knee E > ZEp, D[E) oc E and 
7 = 3.35. 
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In our model the iron spectra are different for the cases 1' {Ec = 3 x 
10*^ GeV), 2' (E,, = 1 X 10^ GeV) and 3' {E^ = 2x 10^ GeV) (see Fig. 4). In 
fact all spectra are not power-law, but in the effective power-law approxima- 
tion they can be roughly characterised by 7 ~ 3.9 for spectrum 1', 7 fa 3.4 for 
spectrum 2' and 7 = 3.3 for spectrum 3'. These spectra, especially 2' and 3' 
are consistent with the Hall diffusion, which predicts 7 = 3.35. 

If galactic sources accelerate particles to energy higher than 1 x 10^ GeV, they 
should be observed in extragalactic background by the "direct" flux. Due to 
multiple scattering of protons in the galactic magnetic field the sources look 
as extensive ones with the typical angular size ~ 20° at distance ~ 10 kpc. 

In conclusion, the model for UHE proton propagation (Section 2) combined 
with measured all-particle spectrum (taken as Akeno-AGASA data) predicts 
the galactic iron spectrum in energy range 1x1 0^ — 1x1 0^ GeV. The predicted 
flux agrees well at ~ 8 x 10"^ GeV and below with the KASCADE data. 
The transition from galactic to extragalactic cosmic rays occurs in the energy 
region of the second knee and is distinctly seen if iron and proton spectra are 
measured separately. In all-particle spectrum this transition is characterised 
by a faint feature because spectral indices of galactic component at < 
1 X 10^ GeV and extragalactic component at > 1 x 10^ GeV are close to 
each other (A7 = 0.25). The spectrum of iron nuclei in 1 x 10® — 1 x 10^ GeV 
energy range agrees with the Hall diffusion. 
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